The astrophysics of planetary systems: formation, structure and 
dynamical evolution 

Proceedings IAU Symposium No. 276, 2010 © 2010 International Astronomical Union 

A. Sozzetti, M.G. Lattanzi & A. P. Boss, eds. DOI: 00.0000/X000000000000000X 

Hamiltonian model of capture into mean 
motion resonance 

Alexander J. Mustill 1 and Mark C. Wyatt 1 

1 Institute of Astronomy, University of Cambridge, 
Madingley Road, CB3 OHA, Cambridge, UK 
email: ajm233@ast.cam.ac.uk, wyatt@ast.cam.ac.uk 

Abstract. Mean motion resonances are a common feature of both our own Solar System and 
of extrasolar planetary systems. Bodies can be trapped in resonance when their orbital semi- 
major axes change, for instance when they migrate through a protoplanetary disc. We use a 
Hamiltonian model to thoroughly investigate the capture behaviour for first and second order 
resonances. Using this method, all resonances of the same order can be described by one equation, 
with applications to specific resonances by appropriate scaling. We focus on the limit where one 
body is a massless test particle and the other a massive planet. We quantify how the the 
probability of capture into a resonance depends on the relative migration rate of the planet 
and particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and 
has decreasing probability for higher eccentricities, although for certain migration rates, capture 
probability peaks at a finite eccentricity. We also calculate libration amplitudes and the offset 
of the libration centres for captured particles, and the change in eccentricity if capture does 
not occur. Libration amplitudes are higher for larger initial eccentricity. The model allows for a 
complete description of a particle's behaviour as it successively encounters several resonances. 
The model is applicable to many scenarios, including (i) Planet migration through gas discs 
trapping other planets or planetesimals in resonances; (ii) Planet migration through a debris disc; 
(iii) Dust migration through PR drag. The Hamiltonian model will allow quick interpretation of 
the resonant properties of extrasolar planets and Kuiper Belt Objects, and will allow synthetic 
images of debris disc structures to be quickly generated, which will be useful for predicting and 
interpreting disc images mad e with ALM A, Darwin/TPF or similar missions. Full details can 
be found in iMustill fc WvattJ (|submittedl) . 
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1. Introduction 

Mean motion resonances (MMRs) occur when two objects' orbital periods are close 
to a ratio of two integers, and a particular combination of orbital angles, the resonant 
argument, is librating. Examples in the Solar System include Neptune and Pluto (3:2 
resonance) and the inner Galilean moons of Jupiter (4:2:1 Laplace resonance). There are 
also now numerous examples of suspected or con firmed MMRs in extrasolar planetary 
systems (e.g., GJ 876 b and c in a 2:1 resonance. iLaughlin fc Chambers! 200 ll ). 

Mean motion resonances also occur between planet s and small dust particles, as seen in 
the Earth's resonant dust ring (jPermott et al.lll994 ). Some extrasolar debris discs, such 
as Ve ga, show evidence of non-axisymmetric clumps ([Holland et al. 1998; IWilner et al 



2002), and several authors have attempted to model these as aris ing from a planet's 
resonant perturbations (e.g., Kuchner fc Holman 2003 ; Wvatt 20031) . 

Resonances are common because convergent migration between orbiting bodies can 
cause them to become captured. There are many mechanisms by which such a semi- 
major axis change can be driven. Early work looked at the tidal evolution of satellite 
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Figure 1. Left: Capture probability for first order resonances, as a function of rescaled ec- 
centricity Jo and migration rate $. Colour scale and contours show capture probability, from 
white (100%) to black (0%). The crosses show the rescaled eccentricity and migration rate cor- 
responding to a particle with Keplerian eccentricity 0.01, migrating at 1 AUMyr" 1 into exterior 
resonances with an Earth-mass planet orbiting a Solar-mass star at 1 AU. The 2:1 and 3:2 
resonances are labelled. Right: Critical migration rate at 1 AU for capture of low-eccentricity 
particles beign trapped into resonances with a migrating planet. At migration rates faster than 
this, capture is impossible at low eccentricity. The critical migration rates for first order reso- 
nances are shown as solid lines. We show typical migration rates for planets embedded in gas 
discs at 1 AU (dashed lines) . Mass-dependent Type I migration occurs for low mass planets and 
mass-independent Type II migration for high mass planets. 



orbits (|GoldreicbJ [r965). In a protoplan etary disc, planets can migrate by tidal interaction 
with the gas disc ( see [C hambers 2009, for a recent review), and small planetesimals by 
aerodynamic drag ( Weidenschillind 1977 ). In a g as-depleted debris disc, planets can mi- 
grate by gravitational scattering of planetesimals ( Kirsh et al.ll2009h . Interpl anetary dust 



drifts towards the Sun under the influence of Poynting- Robertson (PR) drag ([Burns et al 
1979h . 



Resonance capture has been studied by several authors. The regime of adiabatic mi- 
gration, where the migration timescale is much longer than the resonant argument's 
librati on timescale, has been studied extensively a nalytically using a Hamiltonian model 
(e.g., lHenrardl[l98a iBorderies fc Goldreichlll984l) . Rapid migration was studied using 
full N-body models by IWvattl (|2003l henceforth W03) for the case of a plane t migrating 
into a planetesimal disc, and using the Hamiltonian model bv lQuillenl (l2006l henceforth 
Q06) for general migration scenarios. Q06 obtained capture probability as a function of 
migration rate and eccentricity. We have extended this work, and using the Hamiltonian 
model we calculate capture probabilities, libration amplitudes and offsets for particles 
that are captured, and eccentricity jumps for those that pass through the resonance 
without capture. 

The Hamiltonian model we use has several advantages over N-body simulations: (1) 
it allows some results to be derived analytically; (2) it is faster to integrate numerically 
than the 3-body problem; (3) all resonances of the same order reduce to a Hamiltonian 
of the same form, with fewer free parameters than the three-body problem. Once a suite 
of numerical integrations of the Hamiltonian model is performed, it can be applied to 
any system, without the need for running a different N-body integration every time the 
system parameters are changed. 
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2. Model 



We work with a Hamiltonian model of mean motion resonance fe.g-. lMurrav fe Dermott 



1999). We consider the circular restricted three body problem with a massive planet and 



a massless test particle orbiting a central star. The test particle orbits either interior 
or exterior to the planet; the resulting equations of motion are the same, with only a 
change in the scalings between the physical variables and the rescaled variables of the 
Hamiltonian model. The Hamiltonian is 

H = J 2 +ftJ- J k l 2 cos k0, (2.1) 

where 9 is related to the resonant argument and the conjugate momentum J is propor- 
tional to the square of the particle's eccentricity, and k is the order of the resonance. The 
parameter ft measures the proximity in semi-major axis to the resonance, and to simu- 
late migration ft is varied at a constant rate proportional to the physical migration rate. 
The equations of motion arising from this Hamiltonian were integrated numerically. We 
varied the initial momentum Jo and the migration rate ft. There are two parameters in 
the problem: the eccentricity of the particle when it hits the resonance (governed by Jo, 
the initial value of J), and the particle's or planet's migration rate (governed by ft). For 
each point in (Jo, ft) parameter space we integrated 100 trajectories, with the resonant 
argument chosen from a uniform distribution over [0, 2tt). 

Figure 1 (left) shows capture probabilities for first order resonances as a function of 
Jo and ft. We see that capture into resonance is guaranteed for small initial eccentricities 
and migration rates. For low migration rates we are in the well-studied adiabatic regime. 
For low eccentricities capture is certain since the separatrix forms around the initial 
trajectory. For high eccentricities the separatrix forms inside the initial orbit and expands 
to meet it as the migration continues. Capture t hen is probabili stic, with a probability 
that decreases as the initial eccentricity increases (|Henrardlll982l) . For low eccentricities, 
with Jo < 1.3, capture is certain if the migration rate is low and impossible if it is high. In 
the limit of low eccentricity, the transition occurs at a critical migration rate of \ft\ ~ 2.1. 
Certain capture occurs for low migration rates up to Jo ~ 1.3. For higher eccentricities 
(Jo > 1.3), capture is always probabilistic, with a capture probability that is not strongly 
dependent of migration rate; however, if migration rate is too high, then capture is still 
impossible. 

We also found the amplitudes and centres of libration of the resonant argument, and 
the change in eccentricity if a particle is not captured. In the presence of migration, the 
centre of libration is offset from the centre in the absence of migration by an amount 
proportional to the migration rate, almost independent of eccentricity. The amplitude 
of libration increases slightly with migration rate, and significantly with the particle's 
initial eccentricity. If a particle is not captured into resonance, its eccentricity is driven 
down if capture failed d ue to high eccentricity (in agreement with adiabatic theory; 
Murray fc Dermott 19991 ). up if capture failed due to fast migration but low eccentricity, 



and either up or down if both ecce ntricity and migration rate are high. 



In lMustill fc Wvatti ([submitted! ) we also investigated second-order resonances. 
We validated the model against the N-body integrations of W03. The Hamiltonian 
model generally shows good agreement with the results of N-body integrations. 



3. Applications 

Using the Hamiltonian model, the result of an encounter with a resonance in many 
scenarios can be easily determined by rescaling variables and looking up the outcome in 
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the grid of integration results. For example, Figure 1 shows the rescaled eccentricity and 
migration rate for a particle with Keplerian eccentricity 0.01, migrating at 1 AUMyr -1 
into exterior resonances with an Earth-mass planet orbiting a Solar-mass star at 1 AU. 
We see that it is impossible for the particle to be captured into the 2:1 resonance. If 
migration continues it is almost certain, however, to be captured into the 3:2 resonance, 
and if it is not it will be captured into the 4:3 resonance. Because determining the fate 
of a particle is reduced to looking up the outcome in a data table, the fate can be 
determined much faster than it can by performing N-body integrations. This makes the 
model very suitable for generating synthetic images of debris discs with many particles 
(outward migration of the planet and inward migration of the dust can both be handled 
in the same way) , and for determining the outcome of resonant encounters during planet 
formation in population synthesis models. 

3.1. Debris discs 



In iMustill fc Wvattl ([submitted! ) we outline how the degree of dynamical excitation of 
a planetesimal disc can affect the observed structure if a planet migrates into it. More 
dynamically excited discs will show weaker resonant signatur es, since the captu re prob- 



ability is lower and the amplitude of libration higher (see also iReche et al.ll2008f) . Future 
work will explore disc structures in detail using this model. 

3.2. Protoplanetary discs 

In Figure 1 (Right) we show the critical migration rate for capture of particles interior 
to a migrating planet. Also indicated are typical migration rates. It is straightforward to 
determine which resonance the planet will capture particles into. For example, a 10 Earth 
mass planet undergoing Type I migration will capture bodies into the 4:3 resonance, 
while a 1000 Earth mass planet undergoing Type II migration will capture bodies into 
the 5:3 second-order resonance. Future work will address the issue of how accurate the 
Hamiltonian model is at describing the dynamics of two massive planets. 
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